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Abstract 

We consider the constraints of causality and unitarity for the low-energy interactions of protons 
and neutrons. We derive a general theorem that non- vanishing partial- wave mixing cannot be 
reproduced with zero-range interactions without violating causality or unitarity. We define and 
calculate interaction length scales which we call the causal range and the Cauchy-Schwarz range 
for all spin channels up to J = 3. For some channels we find that these length scales are as large as 
5 fm. We investigate the origin of these large lengths and discuss their significance for the choice 
of momentum cutoff scales in effective field theory and universality in many-body Fermi systems. 
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I. INTRODUCTION 



Chiral effective field theory describes the low-energy interactions of protons and neutrons. 
If one neglects electromagnetic effects, the long range behavior of the nuclear interactions 
is determined by pion exchange processes. See Ref. UQ for reviews on chiral effective 
field theory. But there are also systems of interest where momenta smaller than the pion 
mass are relevant. In such cases it is more economical to use pionless effective field theory 
with only local contact interactions involving the nucleons. The pionless formulation is 
theoretically elegant since the theory at leading order is renormalizable and the momentum 
cutoff scale can be arbitrarily large 0-11]. This allows an elegant connection with the 
universal low-energ y p hysics of fermions at large scattering length and other systems such 



as ultracold atoms 
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For local contact interactions the range of the interactions are set by the momentum 
cutoff scale for the effective theory. There are rigorous constraints for strictly finite-range 
interactions set by causality and unitarity. Some violations of unitarity can relax these 
constraints if one works at finite order in perturbation theory or includes unphysical prop- 
agating modes with negative norm. However at some point one must accurately reproduce 
the underlying unitary quantum system by going to sufficiently high order in perturbation 
theory or decoupling the effects of propagating unphysical modes. 

The time evolution of any quantum mechanical system obeys causality and unitarity. 
Causality requires that the cause of an event must occur before any resulting consequences 
are produced, and unitarity requires that the sum of all outcome probabilities equals one. In 
the case of non-relativistic scattering, these constraints mean that the outgoing wave may 
depart only after the incoming wave reaches the scattering object and must preserve the 
normalization of the incoming wave. In this paper we discuss the constraints of causality 
and unitarity for finite range interactions. Specifically we consider neutron-proton scattering 
in all spin channels up to J = 3. 

The constraints of causality for finite-range interactions were first investigated by Wigner 
[ijj]. The time delay between an incoming wave packet and the scattered outgoing wave 
packet is equal to the energy derivative of the elastic phase shift, At = 2hd5/dE. If d5/dE 
is negative, the outgoing wave is produced earlier than that for the non-interacting system. 
However the incoming wave must first arrive in the interacting region before the outgoing 
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wave can be produced. For each partial wave, L, this puts an upper bound on the effective 
range parameter, r^, in the effective range expansion, 

k 2L+1 cot 8 L (k) = -— + -r L k 2 + 0(k 4 ). (1) 

CLl 2 



Phillips and Cohen |15[] derived the causality bound for the S"-wave effective range parameter 
for finite-range interactions in three dimensions. Further work on causality bounds and 
universal relations at low energies have been carried out by Ruiz Arriola and collaborators. 
Constraints on nucleon-nucleon scattering and the chiral two-pion exchange potential was 
considered in Ref . 16] . Correlations between the scattering length and effective range were 



discussed for one boson exchange potentials [17[ and van der Waals potentials 
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The causality bounds for arbitrary dimension d and arbitrary angular momentum L were 



derived in Ref. 
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20(. Let R be the range of the interaction. For the case d = 3, it was 



found that the effective range parameter must satisfy the upper bound 

2T{L-\)Y{L + \) /rx-MH-i 
r L < o L {r) = z — £- y-j 

4 1 /r\2 2vr 1 /r\2£+3 

~T+\1^\2) + r(L + f)r(L + |)4 12J ' (2) 

for any r > R. This inequality can be used to determine a length scale, R b , which we call 
the causal range, 

r L = b L (R b ). (3) 

The physical meaning of R b is that any set of interactions with strictly finite range that 
reproduces the physical scattering data must have a range greater than or equal to R b . 

In this paper we extend the causality bound and causal range and to the case of two 
coupled partial-wave channels. For applications to nucleon-nucleon scattering the relevant 
coupled channels are 3 S\- 3 Di, 3 P2- 3 F 2 , 3 D 3 - 3 G^, etc. As we will show, there is some 
modification of the effective range bound in Eq. (J2J) due to mixing. For total spin J we 
show that the lower partial-wave channel L = J — 1 satisfies the new causality bound, 

2 r(j+i)r(j+|) 



rj-i < &j-i(r) - 2q^ ^ ^ y , (4) 

where go is the first term in the expansion of the mixing angle ej in the Blatt-Biedenharn 



eigenphase convention 2l| . 



tanej(k) = q k 2 + q x k 4 + 0(k 6 ). (5) 



We note that the last term in Eq. (jlj) is negative semi-definite and diverges as r — > 0. From 
this observation we make the general statement that non-vanishing partial-wave mixing is 
inconsistent with zero-range interactions. We will explore in detail the consequences of this 
result as it applies to nuclear effective field theory. 

We also derive a new causality bound associated with the mixing angle itself. Using 
the Cauchy-Schwarz inequality we derive a bound for the parameter q\ in the expansion 
Eq. ([5]). This leads to another minimum interaction length scale, which we call the Cauchy- 
Schwarz range, R . We use the new causality bounds to determine the minimum causal 
and Cauchy-Schwarz ranges for each 2S+1 Lj channel in neutron-proton scattering up to 
J = 3. Since the long range behavior of the nuclear interactions is determined by pion 
exchange processes, one expects R b ~ R c ~ s ~ m^ 1 = 1.5 fm. However in some higher 
partial-wave channels we find that these length scales are as large as 5 fm. We show these 
large ranges are generated by the one-pion exchange tail. Using a potential model we show 
that the causal range and Cauchy-Schwarz range are both significantly reduced when the 
one-pion exchange tail is chopped off at distances beyond 5 fm. We discuss the impact of 
this finding on the choice of momentum cutoff scales in effective field theory. 

In the limit of isospin symmetry our analysis of the isospin triplet channels can also be 
applied to neutron-neutron scattering and therefore has relevance to dilute neutron matter. 
The physics of dilute neutron matter is important for describing the crust of neutron stars 
as well as connections to the universal physics of fermions near the unitarity limit. Our 
analysis of the causality and unitarity bounds show that there are constraints on the universal 
character of neutron-neutron interactions in channels with partial-wave mixing as well as 
higher uncoupled partial-wave channels. In other words some low-energy phenomenology 
cannot be cleanly separated from microscopic details such the range of the interaction. 
Reviews of the theory of ultracold Fermi gases close to the unitarity limit and their numerical 



simulations are given in Re 
length can be found in Ref. 
at unitarity. 
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22]. A gen era 



overview of universality at large scattering 



12l | . See Ref. 23|,|24j for reviews of recent cold atom experiments 
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II. UNCOUPLED CHANNELS 



We analyze in this section the channels with only one partial wave, L. We summarize 
the results obtained in Ref. {20)]. For simplicity, we will assume throughout the calculations 
that the interaction has finite range R, and we use units where H = 1. For the two-body 
system the total wave function in the relative coordinate is 



¥ k \r) = R[ k \r)Y L 



(6) 



Rj/(r) is the radial part, and the Y l,m l {&i 0) are the spherical harmonics. The radial wave 
function R^'(r) satisfies the radial Schrodinger equation, 



1 d 

r 2 dr V dr 



r- 



rR 

Rf\r) + 2fi / W(r,r')U { L k \r')dr' = 0. (7) 
Jo 



Here \x is the reduced mass. We write W(r,r') for the non-local interaction potential as a 
real symmetric integral operator. Since the interaction has finite range R, we require that 
W(r, r') = 0forr>Rorr'>R. U {k) (r) is the rescaled the radial function, 

U?\r) 



R[ k \r) 



kr 



(8) 



The effective range expansion for channels with a single partial wave is, 

k 2L+l cot 5 L (k) = — + l -r L k 2 + O (k 4 ) . 
Ql 2 



(9) 



It is shown in Ref. 



20] that for any r > R the effective range satisfies 



where biir) is 



b L (r) 



tl = b L (r) 



2ir 



pr 


U<£\r>) 
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2 dr' 



(10) 



£r(L + f)r(L + §) \2) 



r \ 2L+3 



1 4 _ ^r_y _ 2T {L-\)T(L + \) ^-2L + i 



a L L + ± \2 



7T 



(11) 



Since the wave function is real and the integral term in Eq. (TTU]) is positive semi-definite, 
this equation puts an upper bound on the effective range, < 6i(r). 
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III. COUPLED CHANNELS 



In this section we derive the general wave functions for spin-triplet scattering with mixing 
between orbital angular momentum L = J — 1 and L — J + 1. The coupled-channel wave 
functions satisfy the following coupled radial Schrodinger equations, 

d 2 



R 



+ 2/i I [W n (r, r^ufl^r') + W 12 (r, r')V}^(r')} dr 
'o 



r(k) 



(12) 



dr 2 



-k 2 + 



(J + l)(J + 2) 



+ 2/i / [W 21 (r, r')Uj\(r') + W 22 {r, r')vj%{r')] dr' = 0. 
Jo 



(13) 



W(r, r') 



(14) 



Here the non-local interaction potentials are represented by a real symmetric 2x2 matrix 
W(r, r'), 

Wn (r, r') W 12 (r,r') 
^W 12 (r,r') W 22 (r,r') 

In Eq. ( fl2|) -( fr3|) the Uf^^r) corresponds with the spin-triplet L = J — 1 channel and the 
Vj+i( r ) is for the spin-triplet L — J + 1. These wave functions are the rescaled form of the 
radial wave functions. In the non-interacting region r > R the coupled radial Schrodinger 
equations reduce to the free radial Schrodinger equations 



d 2 fc2 | MJ ~ I) 
dr 2 r 2 



0. 



cP 



dr 2 



k2+ {J + 1){ j + 2 y ]vj%[u 



0. 



(15) 
(16) 



The solutions of these differential equations are the Riccati-Bessel functions, 

Uf\(r) = A 1 S J -. 1 (kr) + B x Cj^{kr), 

A 2 S J+1 (kr) + B 2 C J+1 (kr). 

where Ai i2 and B\ )2 are amplitudes associated with incoming and outgoing waves, respec- 
tively. More details regarding the Riccati-Bessel functions are given in Appendix IA II The 
relation between incoming and outgoing wave amplitudes is 



(17) 
(18) 



B = KA, 



(19) 



K 1 is the reaction matrix and is defined in terms of the unitary scattering matrix S 

K = i(l-S){l + S)-\ 



(20) 



Therefore, the Eq. ( Tl9|) is written as 



B = SA, 



(21) 



where an d i?i 2 are rescaled amplitudes associated with incoming and outgoing waves. 
For two coupled channels the 2x2 scattering matrix can also be made symmetric. It is 
possible to write several different 2x2 S'-matrices which satisfy the unitarity and symmetry 
properties. In the literature, there are two conventionally used 2x2 S'-matrices j^Ql. In 



this study we adopt the "eigenphase" parameterizations o 



Blatt and Biedernharn 



and 



the relations between the eigenphase and nuclear bar 
Appendix IB1 

The S-matrix can be diagonalized by an orthogonal matrix U 

2iSa Q 



25J parameterizations are shown in 



s d = usu~ l 



e 2iS v 



(22) 



that contains one real parameter e, 



U 



cos e sin e 



— sin £ COS £ 



(23) 



S a (k) and dp(k) are the two phase shifts, and e(k) is the mixing angle. The S-matrix 
explicitly is 

(g2i<5 a g2i<5,3j 



s 



e 2i5 a cog 2 £ _|_ e 2«<5 /3 g j n 2 £ cog £ g - n £ ( e 2i& a 



cos £ sin e e 



^2iS a _ p 2iS g ^ p 2i5, 



sin 2 £ + e 2l& P cos 2 £ 



(24) 



The eigenvalue equation S \X) = A \X) results in eigenvalues Ai = e 2l6a and A2 = e 2l5/3 , 
with corresponding eigenstates, 



l*i> 



COS£ 



sine 



and \X 2 ) 



sine 



COS£ 



(25) 



which satisfy the orthogonality condition 



{X 1 1 X 2 ) = 0. 



(26) 
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We can write Eq. (l2Tj) as 

B2a B 2 /3 

where the matrices A and B are 




(27) 



e lSa cos e — e 4,5,3 sin e 



A = I ' I , (28) 

e a sin £ e 8 cos e 

, e^ Q cose — e i<5fl sine \ 
B =\ • (29) 

e <* sin £ e 3 cos £ / 

We now define some additional notation. We write all a-state phaseshifts 5 a (k) as 
8j-i(k) and all /3-state phaseshifts 5p(k) as The notation is appropriate since in 

the k — >■ limit the a-state is purely L = J — 1 and the /3-state is purely L = J + 1. We 
also drop the superscript k in the wave functions. We choose the normalization of the wave 
function to be well-behaved in the zero-energy limit. Using the relations 

CMkr) — » (frr) JT ' 2 J± T(J± 1 + 1/2), (31) 

as fc— >0 yj 71 

and removing an overall phase factor, we get wave functions of the form 

U a (r)= cos ej(k) k^cot 5j^(k)Sj^(kr) + Cj^(kr)\, (32) 

V a (r)= sine j(k) k^ 1 [cot 5 j^(k)S J+1 (kr) + C J+1 (kr)}, (33) 

U p (r) = -sme.j(k) k J+1 [cot 6 J+1 ( k) Sj-^kr) + Cj^(kr)], (34) 

Vp{r) = cosej(k) k J +1 [cot 5 J+1 (k)S J+1 (kr) + C J+1 (kr)]. (35) 

For later convenience we define 

s J±1 (k,r) = k- J ^- 1 S J±1 (kr), (36) 
cj ±1 (k,r) = k J±l C J±1 (kr). (37) 

Eq. (J36l)- (l371) with Eq. (1Aljl - (lA2jl indicate that s j±i(k, r) and Cj±\(k, r) are analytic functions 
of k 2 , and they can be written as 

sj±i(k, r) = s ,j±i{r) + k 2 s 2 , J±1 (r) + O (k 4 ) , (38) 
cj±i(k, r) = c , J±1 (r) + fc 2 c 2 , J±1 (r) + O (k 4 ) . (39) 
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The explicit form for the terms in these expansions are given in Appendix IA II Therefore 
Eq. ((32D-fl35D become 

U a {r) = cosejik)^- 1 cot5j^{k)sj^{k,r) + cj^{k,r)}, (40) 
V a {r) = smej(k)[k 2J+1 cot5j^(k)s J+1 (k,r) + k- 2 c J+1 (k,r)}, (41) 

U p {r) = -smej(k)[k 2J+1 cot8 J+1 (k)sj^(k,r) + k 2 Cj ^(k,r)}, (42) 
Vp(r) = cose j(k)[k 2J+3 cot 5j +1 (k)s J+1 (k,r) +c J+1 (k,r)\. (43) 

Analogous to the effective range expansion defined in Eq. (jUJ), the two-channel effective 
range expansion has the following form, 

fc£«+§£-ijfe£«+* = L + -r Li .k 2 + C(A; 4 ), (44) 

3.7",.. 2 

where ax y is the scattering length matrix, r^. is the effective range matrix, and k Lij is the 
diagonal momentum matrix diag(A; J_1 , k J+1 ). The two-channel effective range expansion in 
the Blatt and Biedernharn parameterization is 

k Li * + *Trkr x U- x k Li i + * = — + -r L ..k 2 + 0(Jfe 4 ) (45) 

where a^, = diag(aj_i, aj+i), and r£ y = diag(rj_i, O+i)- in addition we get an analytic 



expansion for the tangent mixing angle 26] 



tan ej(jfe) = q k 2 + qi k 4 + O (k 6 ) (46) 

with mixing parameters go an d qi- Now using Eq. (|46p we obtain the following final forms 
of wave functions for r > R, 

U a {r) = — —s 0iJ -i(r) + c a j-i(r) 
aj-i 

+ fc 2 {^O_ lS0 ,j-i(r) - —^^(r) + c 2> j^(r)\ + 0(k 4 ), (47) 
^ a (r) = g c ,j+i(r) + /c 2 {gic 0ij+1 (r) + g c 2 ,j+i(r)} + 0(k 4 ), (48) 

= <?o— ^s .j-i(r) 
a J+i 

+ k 2 {q Q —s 2yJ ^{r) - q JJ±L s ^ x {r) + q x — s 0> j-i(r)} + (9(A; 4 ), (49) 
I a J+ i 2 a J+ i J 



-1 

a J+i 



So,J+i(r) + c 0iJ+ i(r) 



+ k 2 {\r J+l s^ J+l {r) L S2J+1 ( r ) + ca J+1 (r)} + 0(k A ). 

i> Z aj+i J 



(50) 



As in the single channel case, the tool that we use to derive the causality bound is the 
Wronskian identity. Through the derivation we assume that the potential is not singular 
at the origin, and that regular solutions of the Schrodinger equations U(r) and V(r) for two 
different values of momenta, k a and k b , satisfy 



lim U b (p)U' a (p) = lim U a (p)U' b (p) = 0, 

lim .V b (p)V:(p) = lim V a (p)VM = 0. 
For 7 = a, (3 states we obtain 

{kl-kl) / [Ua^U^ + Va^V^r'^dr' 
Jo 

= W[U ai (r), U b ^(r)} + W[V aj (r), H 7 (r)], 
and for the combination of a and (3 states, we get 

(** - k\) [ [U aa {r')U bp {r') + V aa (r')V b p(r') + U ba (r')U a p(r') + V ba (r')V a p(r')] dr' 



(51) 
(52) 



(53) 



= W[U aa (r), U bp {r)\ + W[U a p{r), U ba (r)} + W[V aa {r), V bp {r)\ + W[V af3 {r), V ba (r)}. 

(54) 

The Wronskian of the a-state wave functions and the /3-state wave functions for the non- 
interacting region r > R are given in Appendix IA2I 

In Eq. (|53|) . we set k a = and take the limit k = k b — > 0. In the region r > R we obtain 
the following relations for the effective range parameters, 



rj-i = &j-i(r) + 2g 2 W/[c 2 (r),c (r)] J+1 - 2 f ([C/f (r')] 2 + K°>(r')] 2 ) dr' 



(55) 



O+i = &j + i(r) + 2g 2 ^^[ S2 (r), S0 (r)]j-i-2 / ( feV )1 + k (0) (r')l (56) 



V+i 



Here 6j T i are 



V=Fl 



+ W[s (r), c 2 (r)] JT1 + 2W[c 2 (r), c (r)] Jt1 , 



(57) 
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which reduce to the form 
1 



2tt 



i^ir(j=Fi + |)r(j T i + |) V2 



r x 2(J T l)+3 



1 



■\2 2r(JTl-|)r(JTl + i) /r\-2(^Ti)+i 



7T 



(58) 



In Eq. ( EH) ; we set fc a = and take the same limit, A; = — > 0. In the region r > R we 
obtain 



d.r(r) - 2 



Here dj(r) is 



(59) 



dj(r) = - %- 



-W[s 2 (r), s„(r)] j-i + 2g W/[c 2 (r), c (r)] 



aj-iaj + i 
2 r 

+ go ^ W[c 2 (r), s (r)]j-i - W[c 2 {r), s {r)) J+1 



J+i 



and this can be written as 
-go 



2ir 



aj-iaj+i r (I + j) r (| + j) V2 
go 4 



r \ 2J+1 



r 2 _ 2a r(j+|)r(j+|) ,r V »-i 

a J+1 {2J- l)(2J + 3) yo TT V2, 



(60) 



(61) 



All of equations derived here have been numerically checked using a simple potential 
model. The numerical calculations using delta-function shell potentials with partial-wave 
mixing have been performed, and details are given in Ref. j^ . 



IV. CAUSALITY BOUNDS 

The terms in the integrals in Eq. ( 15 5 j) and Eq. (15 6j) are positive semi-definite since the wave 
functions are real. Therefore Eq. ( 155]) and Eq. (156]) place upper bounds for the effective range 
rj_! and rj +1 respectively. As noted in the introduction, these upper bounds result from 
the causality and unitarity in the quantum scattering problem. Our results are extensions 
of single- channel results in Ref. 15 1 for the S'-wave in three dimensions and in Ref. [l9] for 
arbitrary angular momentum and arbitrary dimensions. 

The causality bounds for the lower and higher partial-wave effective ranges are 

2 r(j + |)r(j+|) 



rj-i < 6j_!(r) - 2q^ ^ ^ [-) , (62) 
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r J+ i < b J+1 (r) + 



71 



t 2 J+1 r(j + l)r(j + 



2J+1 



(63) 



We note that the effective range bounds are modified due to partial-wave mixing. The 
causality upper bound for rj_i is lowered by the negative term on the right hand side of 
Eq. (I62p . while the causality upper bound for the higher partial- wave is increased by the 
term on the right hand side of Eq. f l6"3"|) . When q is nonzero and we take the limit of 
zero range interactions, Eq. (162]) tells us that rj_i is driven to negative infinity for any 
J. We conclude that the physics of partial-wave mixing requires a non-zero range for the 
interactions in order to comply with the constraints of causality and unitarity. In Ref. [ljj] 
a similar negative divergence in the effective range parameter was found for single-channel 
partial waves with L > 0. What is interesting here is that the negative divergence of the 
effective range occurs already in the 3 Si channel due to partial- wave mixing. 

We note that the integral terms in Eq. (|55|) . Eq. (|56|) and Eq. (159|) are closely related. 
Analysis of these equations using the Cauchy-Schwarz inequality provides another useful 
relation for the coupled-channel wave functions. For real functions fi{r), M r )> 9i{ r ) and 
gi (r), the Cauchy-Schwarz inequality is 



\h{r) h(r)} 



Mr) 



dr 



[9i(r) g 2 (r)} 



9i{r) 
92(r) 



dr 



> 



[M r )9i{r) + M r )92{r)] dr 



(64) 



When we apply the inequality to our coupled wave functions, we get 



fj^(r)g J+1 (r) > [hj(r)f 



(65) 



where 



and 



/j_i(r) = 6j_!(r) - 2gj 



,r(j+mj+ 



71 



g J+ i{r) = b J+ i{r) + 



hj(r) = dj(r) - q v 



-2J-1 



2gp 7T /r\ 

« 2 , + ir(J+|)r(j+f) \2J 

2 



2J+1 



rj+i, 



(66) 
(67) 

(68) 



This inequality is used to define a Cauchy-Schwarz range, R c ~ s , as the minimum r for each 
coupled channel where Eq. fl62|) . Eq. fl63l and Eq. fl65|) hold. 
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V. NEUTRON-PROTON SCATTERING 



We now apply our causality bounds to physical neutron-proton data. In this study, we 
use the low energy neutron-proton scattering data (0-350 MeV) from the NN data base by 
the Nijmegen Group 28] . Table |T] - IIHI show the low-energy threshold parameters in the 
eigenphase parameterization for the Nijmll and the Reid93 potentials. These parameters 
are calculated using the results obtained in Ref. |29j for the low-energy threshold parame- 
ters of the nuclear bar parameterization and relations between eigenphase and nuclear bar 
parameterizations given in Appendix [Bj Using these numbers we analyze Eq. fl55|) . Eq. f )56|) 
and Eq. (I5l?j) . as well as causality bounds for the uncoupled channels. 



TABLE I: The eigenphase low energy parameters of uncoupled channels for neutron-proton scattering by 

the Nijmll and the Reid93 interaction potentials. 



Channel 


a L [ftn 2i+1 ] 
Nijmll (Reid93) 


r L [fm- 2L+1 ] 
Nijmll (Reid93) 


So 


-23.727 (-23.735) 


2.670 (2.753) 


l Pl 


2.797 (2.736) 


-6.399 (-6.606) 


3 Po 


-2.468 (-2.469) 


3.914 (3.870) 


3 Pi 


1.529 (1.530) 


-8.580 (-8.556) 


l D 2 


-1.389 (-1.377) 


14.87 (15.04) 


3 D 2 


-7.405 (-7.411) 


2.858 (2.851) 


l F 3 


8.383 (8.365) 


-3.924 (-3.936) 


3 F 3 


2.703 (2.686) 


-9.932 (-9.994) 



A. Uncoupled Channels 

We start with channels of a single uncoupled partial wave. Since there is no mixing 
between different partial-waves, we evaluate Eq. f )55|) and Eq. (1561) with zero mixing angle, 
and we obtain the following equation for the effective range 



r L = b L (r) 



C/(°)(r') 



dr' 



(69) 
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TABLE II: The eigenphase low energy parameters of coupled channels for neutron-proton scattering by 

the Nijmll and the Reid93 interaction potentials. 



Channel 


a L [fm 2i+1 ] 


r L [fm- 2L+1 ] 




Nijmll (Reid93) 


Nijmll (Reid93) 


3 Si 


5.418 (5.422) 


1.7531 (1.7554) 


3 D 1 


6.0043 (5.9539) 


-3.523 (-3.566) 


3 P2 


-0.2844 (-0.2892) 


-11.1465 (-10.7127) 


3 F 2 


8.126 (7.882) 


-5.640 (-5.821) 


3 D 3 


-0.1449 (-0.177) 


288.428 (198.528) 


3 G 3 


648.813 (534.594) 


-0.03306 (-0.0529) 



TABLE III: The eigenphase low energy mixing parameters of coupled channels for neutron-proton 
scattering by the Nijmll and the Reid93 interaction potentials. 



Mixing angle 


qo [fm 2 ] 
Nijmll (Reid93) 


Ql [fm 4 ] 
Nijmll (Reid93) 


£l 


0.303987 (0.303394) 


-2.00228 (-1.99129) 


£2 


-5.65752 (-5.5325) 


65.8602 (64.2979) 


£3 


66.6632 (54.7062) 


340.988 (94.9015) 



where bi is given in Eq. ( fill . These solutions were derived by Hammer and Lee [20j for 
arbitrary dimension and angular momentum. 

Here, we analyze the causality bound of the effective range for L < 3 using the scattering 
parameters in Table [H In Fig. (JT]), we plot — rj\ for all of uncoupled channels with 

L < 3. The physical region corresponds with — > 0. 

For 5-wave scattering 



2 2 
b (r) = 7T-2?" 3 - — r 2 + 2r, 

OQ/ri 



for P-wave, 



bi(r) 



a 



2r 2 



2r £ 



45a 2 3ai 



(70) 



(71) 



14 




FIG. 1: The plot of [i>z,(r) — rz,]/2 as a function of r for neutron-proton scattering via the Nijmll 

potential in the 2S+1 Lj channel. 



for .D-wave, 



for F-wave, 



and for G-wave, 



b 2 (r) 



h(r) 



2 7 2 2 6 

~r r -. 

1575a9 5a2 r d 



2r 9 2r 2 



99225a^ 7a 3 



6 4 (r) 



90 



2r n 2r 2 3150 



9823275a| 9a 4 r 7 



(72) 
(73) 
(74) 



B. Coupled Channels 

We now analyze channels with coupled partial waves. We plot Eq. fl66|) and Eq. ([6] 
for all coupled channels with J < 3. The physical region correspond both /j_i(r) > and 
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FIG. 2: The plot of /j_i(r) as a function of r for neutron-proton scattering via the Nijmll potential for 
J < 3. Here /i(r) is rescaled by a factor of 0.01 and fi{r) is rescaled by a factor of 10~ 4 . 




FIG. 3: The plot of gj+\(r) as a function of r for neutron-proton scattering via the Nijmll potential for 

J < 3. Here gs(r) is rescaled by a factor of 0.1. 
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1. 3 Si- 3 Di Coupling. 



We now consider Eq. (l55j) - (!60|) for the 3 Si- 3 -Di coupled channel. We evaluate the 
Wronskians for J = 1 and get 



bo(r) - ql^ - r = 2 f (^V)]* + KV)]*) dr', 



,2r 3 



W + q>— 2 - r 2 = 2 / ( C/^(r') + Vj>') 



r(0). 



1 2 



'(0), 



1 2 



dr', 



di(r) - gi- = 2 / [^(0^(0 + K to (rO%(0] 

02 



&o( r ) an d b 2 (r) are given in Eq. ( 1701) and in Eq. (172]) . respectively, and di(r) is 

di(r) 



1 2r 3 1 4r 2 6 
-Qo T- + % - 



(75) 

(76) 
(77) 

(78) 



Using the scattering parameters in Table HTl- IHIl we plot Eq. (1751) . Eq. (176j) and Eq. (ITT]) as 
functions of r. In Fig. (j3J) we show the physical region where the causality bounds fo(r) > 0, 
gi(r) > 0, and fo(r)g2(r) > h\(r), are satisfied. Here we have 

6 



fo(r) 
02 (r) = 



JL r 3 _ A r 2 + 2r - ql— - r , 
6clq ao r d 



2 2 6 , 2 2r 3 



1575a 2 



5ac 



r ~ -3 _r 2' 



, . , 1 2r 3 1 4r 2 6 2 

M r ) = 5- + qo — =- - qo^ - qi—. 



(79) 

(80) 
(81) 



2. 3 P 2 - 3 F 2 Coupling. 



In the 3 P2- 3 -p2 coupled channel Eq. (155]) - (1601 take the following forms, 

,90 



h(r)-qf^-n = 2 / ([Ut?\r')] 2 +[ V V\r')] 2 )dr', 



1 2r 5 

W + ^--r 3 = 2 



dr', 



da(r)-gi- = 2 / [L^OMO + V Qa (r')%(r')] dr'. 

a 3 Jo 



(82) 

(83) 
(84) 
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~i 1 1 1 1 1 r 

Re[GMr)g 2 (r)} 

-Re[{f (r)g 2 (r)} m ] 




FIG. 4: We plot Re \/ fo(r)g2(r) , —Re \/ fn(r)g2{r) , and Re[/ii(r)] as functions of r for neutron-proton 

scattering in 3 Si- 3 Di coupled channel. 



h\ (r) and 63 (r) are defined in Eq. flTTj) and Eq. f!73|) . respectively, and ^(r) is 

, . . 1 2r 5 1 4r 2 90 

aia 3 45 a 3 21 r a 

The causality bounds are fi(r) > 0, ^(r) > 0, and fi(r)g 3 (r) > /i 2 (r), where 



AW 



2r 5 2r 2 2 2 90 



03 M 



45aj 3ai r 

2r 9 2r 2 90 2 1 2r 5 
99225a? ~ 7^ ~ + 45~ 



- ^3, 



(85) 



, . , 1 2r 5 1 4r 2 90 2 

/i 2 (r) = -g + 9o— ~ Qo^ ~ qi — ■ 

aia 3 45 a 3 21 r° a 3 

In Fig. (J5J) we show the physical region for the 3 P2- 3 -^2 coupled channel wave functions. 



(86) 

(87) 
(88) 



3. 3 D 3 - 3 G 3 Coupling. 



For J = 3, the 3 D 3 and 3 G 3 channels are coupled. In this case Eq. f l55|) - fl60|) read 
W ~ ^ - ^ = 2 £ {H 0) (r')] 2 + KV)] 2 ) dr\ 



Uf\r') + Vp(r') 



4(r)-gi- = 2 / [U 0a {r')Uw{r')+V Qa {r')V Qp {r')]dr'. 
04 Jo 



(89) 

(90) 
(91) 
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FIG. 5: We plot Re \f fi{r)gz{r) , —Re ^//i(r)g 3 (r) , and Re[/i 2 M] as functions of r for neutron-proton 
scattering in 3 i-2- 3 -F2 coupled channel. The functions are rescaled by a factor of 0.01. 



Here a\(r) is 



d (r) - % 2r? + q 1 ^ g 3150 
a 2 a4 1575 a 4 45 r 7 



The causality bounds are again f2(r) > 0, ^(r) > 0, and f2{r)g±(r) > h\{r), where 



Mr) 



2r 7 2r 2 6 2 3150 
1575a 2 , 5a2 r 3 ^° r 7 



^2, 



04(r) 



2r 



ii 



9823275a? 9a 



2r 2 3150 + q 2 2r 7 



ol 1575 



n, 



q 2r 7 1 4r 2 3150 2 

+ <?o — - % — = — qr 



0,2(14 1575 a 4 45 



ai 



(92) 

(93) 

(94) 
(95) 



We show plots for the 3 Ds- 3 G^ channel in Fig. (|6]). 
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Re[{/2(^W) n ] 
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FIG. 6: We plot Re V^Wff4W > ~^- e V h( r )94{ r ) , and Re[/i3(r)] as functions of r for neutron-proton 
scattering in 3 Z?3- 3 G3 coupled channel. The functions are rescaled by a factor of 0.01. 



VI. RESULTS AND DISCUSSION 

In this section we present the results for the causal and Cauchy-Schwarz ranges, R b , and 
R c ~ s . We use the Nijmll scattering data for neutron-proton scattering presented above. In 
Table IIVI we show results for the causal range for all uncoupled channels by setting 



br,(r). 



(96) 



In Table IVl we determine the causal range for all coupled channels using Eq. ( |66l) -(l671). Also, 
we find the Cauchy-Schwarz ranges shown in Table |Vl] using Eq. ( 1651) . We find that in some 



TABLE IV: The causal ranges for uncoupled channels. 



Channels 


'So 


'Pi 


3 p 

Pq 


3 Pi 


l D 2 


3 D 2 


^3 


3 F 3 


R b [fm] 


1.27 


0.31 


3.07 


0.23 


3.98 


4.91 


1.88 


1.56 



TABLE V: The causal ranges for coupled channels. 



Channels 


3 5x 


3 Dt 


3 P 2 


3 F 2 


3 D 3 


3 c 3 


R b [fm] 


1.29 


1.20 


2.23 


1.73 


4.03 


3.92 
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TABLE VI: The Cauchy-Schwarz ranges for coupled channels. 



Channels 


3 S 1 - 3 D l 


3 P 2 - 3 F 2 


3 D 3 - 3 G 3 


R c ~ s [fin] 


1.29 


4.65 


5.68 



channels the causal and Cauchy-Schwarz ranges are surprisingly large, and it is worthwhile 
to probe the origin of these large ranges. 

It is convenient to collect together some of the key formulas derived above. The Cauchy- 
Schwarz inequality has the form 



fj-i(r)gj + i(r)>[hj(r)} 2 , 



(97) 



where 



f M h M 9 2^+1)^+1) -2.7-1 



7T 



-O-i, 



= 6j+i(r) + 



7T 



r \ 2J+1 



a 2 J+1 r(j+|)r(j+|) 

2 

M r ) = - Qi , 

a-j+i 



-O+i, 



(98) 
(99) 

(100) 



Am 



2tt 



(£) 



2(J=Fl)+3 



a2 Tl r(jT-i + f)r(JT-i + f) 

2 2r (j T i-i)r(j^i + i) /r\ fj 

2, 



1 4 / r \ 2 
~ a jTl J=f 1 + | V2/ 



7T 



;ioi) 



dj(r) 



-go 



2tt 



r x 2J+1 



aj_iaj + i r (i + J) r (f + J) 

go 4 



(5) 



+ 



a J+1 (2J-l)(2J + 3) 



2 n r(j + |)r(j + f) 

r 2 - 2g — — — 

7T 



-2J-1 



(102) 



We note that the leading power of r in g J+1 (r) is 

1 2tt /r \2J+5 



(2) ■ (103 > 



4+i r (J + 1 + 1) r ( j + 1 + 1) 

This has a very small numerical prefactor multiplying aJ 2 1 r 2J+5 . For J = 1 the factor 
is 2/1575, for J = 2 it is 2/99225, and for J = 3 it is 2/9823275. Therefore the term is 
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negligible unless r is large compared with (aj +1 ) 1 /( 2J+3 \ If we neglect this term, then the 
term with the leading power of r on the left hand side of Eq. (1971) is the same as that on the 
right hand side, 



2tt 



a^rfj-i + Drfj-i + f) V2 



71 



a 2 J+1 r(j + i)r(j + |) V2; 



r \ 2J+1 



-go 



2vr 



aj_ia J+1 r (| + J) T (| + J) V2 



r\ 2J+1 



(104) 



As a result the curves for fj_i(r)gj + i(r) and [hj(r)] 2 are approximately parallel for large r 
until the term that we have neglected becomes significant. These nearly parallel trajectories 
inflate the value of the Cauchy-Schwarz range r = R c ~ s where the two curves cross. 

For the 3 Si- 3 Di coupled channel we find R is about the same size as the Compton 
wavelength of the pion, m" 1 = 1.5 fm. This is also comparable to what one expects for 
the range of the nucleon-nucleon interaction. However the results are more interesting for 
J > 1. In the 3 P 2 - 3 F 2 channel we have R = 4.65 fm. And for the 3 D 3 - 3 G3 coupled 
channel we find R c ~ s = 5.68 fm. These values are surprisingly large in comparison with 



VII. ONE-PION EXCHANGE POTENTIAL 

We note that there are some channels where the causal range R b is also quite large. By 
definition R c ~ s > R b and so the Cauchy-Schwarz range will then also be large. The causal 
range is the minimum value for r such that 

2 r(j+i)r(j + |) fr^J-i 



/j_i(r) = 6j_i(r) - 2q^ ^ ^ y - rj _i > (105) 

for the lower partial wave, or 

JJ+l(r) = Wr) + g- r(J + ,; r(J + |) (0 2J+1 - r J+1 > (106) 

for the higher partial wave. For uncoupled channels we take go — 0. 

The largest values for R b occur when the effective range parameter is positive or near 
zero. See for example the causal ranges for the 1 D 2 , 3 D 2) 3 D 3 , and 3 G^ channels. What 
happens is that the function /j_i(r) or gj + i(r) remains negative with a rather small slope 
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until r becomes quite large. The small slope is again associated with the fact that the term 
with the highest power of r has a small numerical prefactor. 

The range of the interaction plays the dominant role in setting the causal range. In the 
language of local potentials, this is the radius at which the magnitude of the potential is 
numerically very small. However there is also some influence of the exponential tail of the 
potential upon the causal range. 

In all channels where the causal range is unusually large, X D 2 , 3 D 2 , 3 D^, and 3 G%, we find 
that the tail of the one-pion exchange potential is attractive. At smaller radii, the potential 
crosses over at some classical turning point to become repulsive. See for example Fig. 2-4 
in Ref. [30]. 



The detailed mechanism requires further study, but it appears that this geometry can 
cause a near-threshold wavepacket to reflect before reaching the classical turning point, thus 
mimicking a longer range potential. However some fine tuning is needed to produce a large 
causal range, as there is no enhancement in the 1 5'o and 3 Si channels and a smaller amount 
of enhancement in the 3 Po channel. 

There seems to be no such enhancement of the causal range in the Pi, 3 Pi, and 3 D\ 
channels where the tail of the potential is repulsive. In fact, the causal range for the 1 P\ 
and 3 Pi channels are unusually small. This appears be related to quantum tunneling into 
the inner region where the potential is attractive. 

In the following analysis we will investigate the importance of the tail of the one-pion 
exchange potential plays in setting the causal range, R b . We show that even though the 
one-pion exchange potential is numerically small at distances larger than 5 fm, chopping off 
the one-pion exchange tail at such distances produces a non-negligible effect. The one-pion 
exchange potential tail appears to be the source of the large values for R b in higher partial 
waves where the central one-pion exchange tail is attractive. 

If we neglect electromagnetic effects, then the neutron-proton interaction potential at 
long distances is governed by the one-pion exchange (OPE) potential, which in configuration 
space is 

V OPE (r) = V c (r) + S 12 V T (r). (107) 
Here Vc{r) is the central potential, 
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FIG. 7: Plot of the potential matrix elements W n (r) = V c (r), Wi 2 (r) = W 2 i(r) = y/SV T (r) and 
W / 22( r ) = Vc( r ) — 2Vr(r) as a function of r in the 3 Si- 3 -Di coupled channel. 



Vj>(r) is the tensor potential, 

W = fg (^)' <* ■ *) (l + ^ + 5^)^. (108) 
and is the tensor operator, 

5*12 = 3((fl • r)((j 2 ■ f) - (7i ■ (72. (110) 

Here m n is the pion mass, is the nucleon mass, and g n N = 13.0 is the pion-nucleon 
coupling constant. The one-pion exchange potential is local in space and the interaction 
matrix in Eq. (1141) takes the following form for J = 1, 

wf * ( Vc{r) ^ T(r) \x( * 

W(r,r)= \5{r-r). (Ill) 

^ VSV T (r) V c (r) - 2V T (r) J 

In Fig. © we plot W u (r) = V c (r), W 12 (r) = W 2l (r) = y/SV T (r), and W 22 (r) = V c (r) - 
2Vt(t) in the 3 S\- 3 Di coupled channel. 

To demonstrate the origin of large causal ranges found in Table [IV] and Table |V] we 
will present some simple but illustrative numerical examples. For each channel we add a 
short range potential to the one-pion exchange potential in order to reproduce the physical 
low-energy scattering parameters. The specific model we use for the short range potential 
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is not important to our general analysis nor is it the most economical. We choose a simple 
scheme which consists of three well-defined functions in three different regions and which is 
continuously different iable everywhere. The potential has the form 



V(r) = V Gauss (r)9(R Gauss - r) + V Sp i inc (r)9(r - i?G auss )0(i2Exch. - r) 
+ V Exc hXr)9(r - -R Ex ch.), 

where 6 is a unit step function. 

The short-range part is a Gaussian function 



V Gauss (r) = C G e- m o r \ 



112) 



(113) 



The intermediate-range part of the potential is a cubic spline use to connect the short- and 
long-range regions, 

^spiine(r) = C X + C 2 r + C 3 r 2 + C 4 r 3 . (114) 

The long-range part consists of the usual one-pion exchange potential together with two 
additional heavy meson exchange terms, 



W(r) = V£ A ' B (r) + S 12 V^ F {r). 



The central part of the potential is composed of Yukawa functions 



(115) 



9l N ( m n \ 2 J r e 



-mgr 



12tt \ 2M N J | ; 
and the tensor part of the potential has the form 

3 3 



+ C, 



+ c 



B~ 



(116) 



2 / \ 2 



12tt V 2M N 
3 



1 + 



m n r (m^r) 2 r 



D 



mr,r (m^r)' 



1 + 



rriFr 



+ 



(mpr) 



-mpr 



117) 



Here C w = (2S 2 - 3)(2T 2 - 3), C n = 2T 2 - 3, g wN = 13.0, = 140 MeV, and M N = 938.0 
MeV. The coefficients not part of the one-pion exchange potential are used as free parameters 
to reproduce the physical low-energy scattering parameters. Due to the abundance of free 
parameters, the fit process is not unique. But in each case we attempt to qualitatively 
reproduce the shape of the Nijmegenll potentials [30j. And in each case the heavy meson 
masses are kept significantly larger than the pion mass. In Fig. ([H|) we show the potential 
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FIG. 8: Plot of the model potential in the 1 P 1 channel. In this channel S = 0, T = 0, CV = 9, C A = 405.9, 
Cb = 769.5, Cn = 0, Co = 0, Cp = 0, Cq = —7.6, = 10.0m,,., mg = 5.45m w , mg = 4.46771^, 
Ci = -7.464, C 2 = 0.179, C* 3 = 73.933 and C 4 = -78.011. We use R Gauss = 0.2 fm and i? Exch . = 0.7 fm. 

in the 1 P\ channel. Fig. shows the potential in the l D 2 channel, and Fig. fllOp shows 
the potential in the 3 Z?2 channel. 

After having recovered the physical low-energy scattering parameters, we now multiply 
an additional step function to the potential, 



which removes the tail of the potential beyond range R. We then recalculate the low-energy 
scattering parameters with this modification. The results are shown in Table IVHI for the 
1 Pi, 1 D 2 , and 3 D 2 channels. The causal ranges for the 1 D 2 and 3 D 2 channels are quite large 
for the physical scattering data, 4.0 and 4.9 fm respectively. But if we remove the tail of 
the model potential at R = 5 fm, the causal ranges drop to 2.4 and 2.7 fm respectively. 

The tail of the model potential is dominated by the one-pion exchange potential. Even 
though the numerical size of the one-pion exchange potential is small at distances of 5 fm, 
these numerical results show clearly that the one-pion exchange tail is controlling the size of 
the causal range. The one-pion exchange potential tail appears to be the source of the large 
values for R b in higher partial waves where the central one-pion exchange tail is attractive. 



V(r) -»■ V(r)0(R-r), 



(118) 
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FIG. 9: Plot of the model potential in the 1 D 2 channel. In this channel S = 0, T = 1, C„ = —3, 
Ca = —60.5, Cb = —30.0, Cn = 0, Co = 0, C F = 0, Cq = 1.01, m A = Qm^, mg — 6m^, wiq — 7.02m w . 
Ci = 1.463, C 2 = -7.332, C 3 = 10.384 and C 4 = -4.585. We use i? Gauss = 0.2 fm and i? Exch . = 0.9 fm. 
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FIG. 10: Plot of the model potential in the 3 D 2 channel. In this channel S = 1, T = 0, = -3, 

C A = -660, C B = -1140, C n = -3, C D = -930, C F = -927, C G = 2.53, m A = m D = 8™^, 
m B =m F = 5m^, m G = 3.25m„, C x = 2.979, C 2 = -4.026, C 3 = -2.455 and C" 4 = 2.409. We use 

-^Gauss = 0.3 fm and i? Exc h. = 1-4 fm. 
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TABLE VII: The potential range dependence of the causal range in various channels. 



Potential range 

Causal range 


R = 2im 


R = 5 fm 


R = 12 fm 


R = 15 fm 


R = 50 fm 


R\ Pl 


0.4 


0.4 


0.3 


0.3 


0.3 




2.0 


2.4 


3.8 


4.0 


4.0 


R *D 2 


1.3 


2.7 


4.7 


4.9 


4.9 



VIII. SUMMARY AND CONCLUSIONS 



In this paper we have derived the constraints of causality and unitarity for neutron-proton 
scattering for all spin channels up to J = 3. We have defined and calculated interaction 
length scales which we call the causal range, R b , and the Cauchy-Schwarz range, R°~ s . The 
causal range is the minimum value for r such that the causal bounds, 

/,-, w = *,-,m - ^ J+ iw+i) _ rM > 0> (119) 

9j+li r ) = bj+l{ r) + -^- r(j+mj+i) (-) - rj+1 >0, (120) 
are satisfied. For uncoupled channels these bounds simplify to the form 

h(r) = g L (r) = b L (r)>0. (121) 

For coupled channels the Cauchy-Schwarz range is the minimum value for r satisfying the 
causal bounds as well as the Cauchy-Schwarz inequality, 

fj-i(r)g J+ i(r) > [hj(r)} 2 . (122) 

If one reproduces the physical scattering data using strictly finite range interactions, then 
the range of these interactions must be larger than R b and R c ~ s . From these bounds we 
have derived the general result that non- vanishing partial-wave mixing cannot be reproduced 
with zero-range interactions. As the range of the interaction goes to the zero, the effective 
range for the lower partial-wave channel is driven to negative infinity. 

This finding has consequences for pionless effective theory where the range of the inter- 
actions is set entirely by the value of the cutoff momentum. If the cutoff momentum is 
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too high, then it is impossible to obtain the correct threshold physics in coupled channels 
without violating causality or unitarity. In some channels we find that the causal range and 
Cauchy-Schwarz range are as large 5 fm. We have shown that these large values are driven 
by the tail of the one-pion exchange potential. In these channels the problems will be even 
more severe, and the cutoff momentum will need to be rather low in order to reproduce the 
physical scattering data in pionless effective field theory. How low this cutoff momentum 
must be depends on the particular regularization scheme. 

We should note that all of these mixing observables are non-vanishing only when one 
reaches higher orders in the power counting expansion, and there is no direct impact on 
pionless effective field theory calculations at lower orders. See, for example, Ref. js] for 
details on power counting in pionless effective field theory. In the zero-range limit, the term 
which drives the negative divergence of the effective range parameter rj_i is 

nj+iw+l) (iy»-\ (123) 



71 \2, 

At leading order there is no divergence since there is no partial wave mixing and qo = 0. If 
higher-order terms are iterated non-perturbatively as in Ref. [||, then the divergence appears 
at order Q 2 , the first order at which q is non- vanishing. If higher-order terms are iterated 
order- by-order in perturbation theory, then the term in Eq. ( I123[) appears at order Q 4 . This 
is one order higher than the analysis presented in Ref. {9J, and we predict that zero-range 
divergences in rj_\ will first appear at this order. 

It important to note that if one works order-by-order in perturbation theory, then the 
constraints of causality and unitarity always appear somewhat hidden. At every order in 
the effective field theory calculation there are new operator coefficients which appear and 
are determined by matching to physical data. There are no obstructions to setting these 
operator coefficients to reproduce physical values. 

It is only when one iterates the new interactions, i.e., by solving the Schrodinger equation, 
that non-linear dependencies on the operator coefficients appear. In this case one finds that 
the constraints of causality and unitarity give necessary conditions for keeping the operator 
coefficients real. Once we fix the regularization, the bound corresponds with branch cuts of 
the effective theory when viewed as a function of physical scattering parameters. 

These branch cuts cannot be seen at any finite order in perturbation theory. However 
a nearby branch point may spoil the convergence of the perturbative expansion. In this 
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context, our causality and unitarity bounds can be viewed as setting physical constraints 
for the convergence of perturbative calculations in pionless effective field theory. 

If the cutoff is taken too high, a branch cut develops which jeopardizes the convergence 
of the perturbative calculation. Similarly if one does calculations using dimensional regular- 
ization, then the renormalization scale sets the scale at which the infrared and ultraviolet 
Physlcs are related 0. SMar p robl e ms with pe rturba ,ve co„ce would a rls e if 
the renormalization scale is taken too high. 

There is much theoretical interest in the connection between dilute neutron matter and 
the universal physics of fermions in the unitarity limit 
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38j. In the limit of isospin 



symmetry our analysis of the isospin triplet channels can be applied to neutron-neutron 
scattering in dilute neutron matter. In this paper we have shown there are intrinsic length 
scales associated with the causal range and the Cauchy-Schwarz range. When the average 
separation between neutrons is smaller than these length scales, one expects non-universal 
behavior controlled by the details of the neutron-neutron interactions. For the 1 S'o channel, 
R b = 1.3 fm. For the 3 P 2 channel, R b = 2.2 fm, and for the 3 F 2 channel, R b = 1.7 fm. 
For 3 P2- 3 -p2 mixing, we find R c ~ s = 4.7 fm. We see that the physics of 3 P2- 3 -^2 mixing 
will become non-universal at lower densities than the 1 S'o interactions. In particular the 
densities where 3 P 2 superfluidity is expected to occur will be well beyond this universal 
regime. 
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Appendix A: Functions 



1. The Riccati-Bessel Functions 

Sj±i(r) and Cj±\(r) are Riccati-Bessel functions of the first and second kind, which are 
defined in terms of the Bessel functions as 

Sj±i(r) = J—Jj± 1+ i(r) = V^r J±1+1 Y — —=- ~-r 2n , (Al) 

j±i\ ) y 2 j±i+^v i v n+1 r(n + J±l + | v ; 

n=0 \ / \ z > 



0F V 2 } 1 2 ; ^I> + 1)I>-J + 1 + 1/2) 1 ; 

The s(r) and c(r) functions in Eq. ( I38jl and Eq. ( 139]) are written from Eq. (lAlj) and 
Eq. (JX2J), 

,/tt / r \J±l + l 

°°.^) = r(J I 1 + |) ( 5 ) • (A3) 



s 2J±1 r = —^1 (') (A4) 

' v ; r(j±i + |) \2J 

:1 + §) (r 
^ \2 

(r 

^ \2 



WrH EE±I±i) f n--\ (A5) 



T (J ±1- -) /rW=Fi+2 
Wr) = 1 ^ 2) C-) ■ (A6) 



2. Wronskians of Wave Functions 

Here we calculate Wronskians of the U(r) and V(r) wave functions, and Wronskians of 
all possible combination of the so(r), S2{r), Co(r) and C2(r) functions. Wronskians of U a (r) 
and V a (r) for the non- interacting region r > R are 

W[C/ aQ (r),C/ 6a (r)] = (A£ - fe 6 2 )|irj_iW[so(r),co(r)]j_i 



2 

+ — ^h(r), s (r)]j_i + —W[co{r), s 2 (0] j-i 
a j-x a J-i 

+ J— W [s (r),c 2 (r)}j^ + W[c 2 (r),c (r)]jA 
a j -i } 

+ 0(ki)+0(kt), (A7) 
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W[V aa (r), V ba (r)} = (k 2 a - k 2 b )q 2 W[c 2 (r), c (r)] J+1 + 0(k 4 a ) + 0(k A b ). (A8) 
Wronskians of the /3-state wave functions are 



1 



W[U a p(r), U bp {r)} = (k 2 a - k 2 b )q 2 —W[s 2 (r), s „(r)] j-i + O(fc^) + 0(*; 4 ), (A9) 

"j+i 



W[^»(r),^(r)] = (k 2 a - k 2 ){K J+1 W[s (r),c (r)] J+1 



+ -^—W[s 2 (r), s {r)] J+1 + W[co(r), s 2 (r)] J+1 

a J+1 aj+i 



—W[s (r), c 2 (r)] J+1 + W[c 2 (r), c (r)] J+1 ) 

+ (MO) 



+ 

Oj+l 



Wronskian of the combinations of the a and /3-states are 



W[U aa (r), U b p{r)] = q —W[c (r), s (r)} j_ x - k 2 a {q W[s 2 (r), s (r)} 3 _ x 

- q—W[c 2 (r), s {r)]jA + k 2 Jq W[s 2 {r), s (r)]j_i 

- q—W[s 2 (r), co(r)] j_i + g^W[s (r), co(r)] j_i 
-gi^^[so(r),c (r)] J _ 1 ) + 0(A; 4 ), (All) 



^[^(r), C/ 6a (r)] = -q—W[c (r), s (r)\ j_ x + k 2 \q W[s 2 (r), s (r)] j_i 

- q-^— W[c 2 (r), s (r)]j_i j - fc*{g W[s 2 (r), s (r)]j-i 



- 9— W[s 2 (r), c (r)] j_! + g^W%(r), c (r)] ^ 
aj+i 2 



-gi^^[so(r),c (r)] J _ 1 ) + 0(A; 4 ), (A12) 



^[Ka(r),H^(r)] = - g J_^[co(r),s (r)] J+1 W[cb(r),s (r)] 

-W[c2(r), s (r)]j+i - qW[c 2 (r), c (r)] J+1 ) 

L J 



+ k 2 Jq r -^-W[c (r),s (r)} J+1 - q—Wic (r), s 2 (r)] J+1 
12 a J+ i 

+ g^[c (r),c 2 (r)] J+1 } + G(A; 4 ), (A13) 
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W\V aP (r), V ba (r)\ = q—Wic (r), s (r)] J+1 + k 2 J qi — W[c (r), s Q (r)\ J+1 

a .j+i L a .j+i 

+ g— !— W[c 2 (r), s (r)] J+1 - qW[c 2 (r),c (r)}j + i} 

~ k 2 a \q r -^W[c (r),s (r)] J+1 - q —W{c (r), s 2 {r)} J+1 
i l aj+i 

+ qW[c (r),c 2 (r)} J+1 } + 0(k 4 ). (A14) 

We now calculate Wronskians of the so(r), s 2 (r), c (r) and c 2 {r) functions. We find 

^Mr). aa (r)] J = - r(j + J)r(§ + J) ( 5 ) , (A15) 

H'[»o(r),c b (r)] J = -l, (A16) 

2 

W[ao(r) ) c i! (r)] J = -^^ ) (A17) 

2 

^(rJ.cbCr)]^^^, (A18) 

4 

^■^' WTiHi 1 (A19) 

W[q,(r), C2(r)]j = _ £(zi±4£(i±Z) (!) . (A20 ) 
It should be noted that W[/(r), #(r)] = -W[/(r),^(r)]. 



Appendix B: Relations between Eigenphase and Nuclear Bar Parameterizations 



The scattering matrix in terms of the eigenphase parameters was given in Eq. (1211) . The 
scattering matrix in terms of the nuclear bar parameters is 

e 2i ^cos2e ie*^ 8 ") sin 2e , 
S= \ f - - \ \ (Bl) 

ie l( < Sa+s e) sin 2e e 2iS ? cos2£ 



Here 5 a , 8p and e are the nuclear bar phase shifts and mixing angle 
between the eigenphase and the nuclear bar parameters are 



25| . The relations 



sin(5 a - 6,) = (B2) 
sin ze 

S a + S^ = 5 a + S , (B3) 

tan2e , 

tan2e = = =—. (B4) 

sin(5 a - 5 P ) 
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The two-channel effective range expansion is defined slightly differently in the eigenphase 
and the nuclear bar parameterizations. In the eigenphase parameterization, 



k L v+*UK- 1 U- 1 k Li * + * = ~ + \r l3 h 2 + 0{k 4 ), (B5) 
and in the nuclear bar parameterization, 



CLij 2 



k^HK' 1 ^^ = -— + -f t ,k 2 + O(^). (B6) 

otj 2 

Therefore, by straightforward calculations we find the following relations among the thresh- 
old scattering parameters, 

a a =a a , (B7) 

r „ =,„ + m + (B8) 

ap =ap - (B9) 
rp =fp, (BIO) 
Qo =-, (Bll) 

?i = ^ • (B12) 

For the uncoupled channels go and q± are zero, and these relations become a a = a a , r a = f a , 

ap = ap, and rp = fp. 
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